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WAVEFRONT METHOD FOR DESIGNING OPTICAL ELEMENTS 

BACKGROUND OF THE INVENTION 

The present invention relates to the design of optica! elements, in general, 
and to a method for designing optical elements using wavefront calculations, in 
particular. 

Most of the designs for progressive (multifocal) ophthalmic lenses are 
based mainly on the geometrical shape of one or more of the lens surfaces, as 
described, for example in U.S. Patent No. 3,687,528 to Maitenaz. The geometrical 
properties, however, are only indirectly related to the lens 1 actual optical 
performance. For example, the curvature of any of the lens surfaces is only 
approximately related to the lens optical power. In some cases, for example U.S. 
patent Application Serial No. 09/262,341 to Katzman et a!., entitled "METHOD 
FOR THE DESIGN OF MULTIFOCAL OPTICAL ELEMENTS", filed March 4, 
1999, the design is examined for its actual optical performance, and parameters of 
the design process can be adjusted in response to the optical performance, but 
this makes the design process indirect and time-consuming. 

Another typical approach in lens design is to use software means to trace 
a large number of rays through the lens and to deduce from the incident and 
refracted rays the classical aberrations of the lens. The surfaces of the lens are 
represented in terms of a fixed number of basic surfaces such as spheres, 
aspherics, tones and cylindrical surfaces. The ray tracing of a large number of 
rays requires a considerable computational effort, and the limited family of 
surfaces in which one optimizes the lens performance greatly restricts the 
optimization capabilities. 

Several methods were proposed to overcome the undesirable constraint 
of a limited family of optical surfaces. U.S. Patent No. 4,613,217 to Fueter et al. 
proposes to represent an optical surface by splines. U.S. Patent No. 5,886,766 to 
Kaga et al. similarly proposes a progressive ophthalmic lens consisting of three 
portions, at least one of which is divided into smaller pieces that are connected 
together by requiring that the surface is at least twice continuously differentiate 



1 



WO 01/92948 



PCT/IL01/00492 



along the interfaces joining two pieces. The spline representation, while offering 
an advantage over the limited class of surfaces mentioned above, has several 
drawbacks: the strict smoothness requirements along boundaries limits the space 
of parameters for the optimization process, the restriction to rectangular pieces, 

5 and thus to rectangular domains, limits the flexibility of the design, and it is not fully 
natural to prescribe boundary conditions on the lens surfaces and their slopes. 
U.S. Patent Application Serial No. 09/262,341 to Katzman et al. describes a finite 
element method for the surface representation. The finite elements are patches of 
arbitrary polygonal shape, thus yielding flexible designs. The method of U.S. 

10 Patent Application Serial No. 09/262,341 to Katzman et al. requires less 
differentiability along the lines joining the patches than that of U.S. Patent No. 
5,886,766 to Kaga et al. 

The common practice in the ophthalmic industry is to design and 
manufacture semi-finished lenses. The lens is manufactured with one surface that 

15 is given, for example a spherical or toric surface, and one surface that is designed. 
The given surface is later processed to meet the specific prescription of the client. 
Consequently most designs are limited in that they take into account only one 
surface of the lens, assuming, in general, a given spherical or toric other surface. 
Furthermore, the capability of improving the lens performance by designing both 

20 surfaces is not utilized. Recent technological developments enable better control 
of the manufacturing of both surfaces of ophthalmic lenses. Indeed, U.S. Patent 
No. 5,771,089 to Barth, U.S. Patent No. 5,784,144 to Kelch and U.S. Patent 
Application Serial No. 09/262,341 to Katzman et al. propose designs in which both 
the front and back surfaces of the lens have flexible shapes. Further technological 

25 developments enable manufacturing of multi-surface ophthalmic lenses. U.S. 
Patent Application Serial No. 09/262,341 to Katzman et al. proposes designs in 
which several of the surfaces have flexible shapes. 

In the design of multifocal ophthalmic lens, one usually selects the front 
(far from eye) surface to be progressive, while the back (close to eye) surface is 

30 either spherical or toric, where a toric surface might be needed to correct 
astigmatism. U.S. Patent No. 2,878,721 to Kanolt and U.S. Patent No. 6,019,470 
to Mukaiyama et al. disclose ophthalmic lenses in which the back surface is a 
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composition of a progressive surface and a tone surface. This composite design 
has a drawback that each of the two surfaces (toric and progressive) is designed 
separately, and thus the optical behavior of one of the surfaces may conflict with 
the optical behavior of the other. It is therefore beneficial to develop a method for 
5 designing a lens consisting of integral surfaces rather than composite surfaces. 
Such a method could then simultaneously consider all desired optical behavior. 

The actual performance of an ophthalmic lens depends not just on the 
lens itself, but on the full eye-plus-lens system. This becomes particularly 
important, for example, when the lens user suffers from astigmatism and/or 

10 presbyopia. Astigmatism is a condition in which the eye focuses differently in 
different directions. Presbyopia is a condition in which the eye loses some of its 
ability to accommodate, i.e. to focus sharply at nearby objects. The cu rvature of 
the lens of the eye changes as the eye focuses on objects at different distances 
from the eye. As people age, their eyes become less elastic and therefore can 

is change the curvature of the lens only to a certain degree. 

The article by J. Loos, Ph. Slusallek and H.-P. Seidel, entitled "Using 
Wavefront Tracing for the Visualization and Optimization of Progressive Lenses", 
Computer Graphics Forum , vol. 17, no. 3, 1998, pp. 255-264, discusses the 
possibility of accounting for the eye structure in the design process of progressive 

20 lenses. However, Loos et al. do not require a full analysis of astigmatism and 
presbyopia since they assume a spherical or toroidal back surface and optimize 
only the front surface of the progressive lens. 

An important characteristic of a lens is its prism. Prism measures the 
change in the direction of light rays as they are refracted by the lens. Essentially 

25 every lens has some level of prism. Sometimes a lens is processed in order to 
induce some desired prism for a variety of purposes. However, design methods 
taking prism into account are not known. 
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SUMMARY OF THE INVENTION 

There is provided in accordance with an embodiment of the present 
invention a method for designing at least one surface of an optical element. For 
each particular surface, a representation is created with a set of discrete points. 
5 Then for each discrete point in the set, nearest neighbor points from the set are 
selected, the particular surface is approximated in a vicinity of the discrete point by 
a polynomial of a predetermined order in two variables of the particular surface, 
and the coefficients of the polynomial are determined according to the selected 
nearest neighbor points without requiring continuity between polynomials for 
10 neighboring discrete points. 

There is also provided in accordance with an embodiment of the present 
invention a method for designing at least one surface of an optical element. The 
method includes the steps of representing the at least one surface using the 
method of unconstrained patches, choosing a function in parameters of the at least 
15 one surface, and optimizing the function with respect to the parameters. 

There is also provided in accordance with an embodiment of the present 
invention a method for designing at least one surface of an ophthalmic lens. The 
method includes the steps of representing the at least one surface with 
parameters, choosing a function in the parameters, and optimizing the function 
20 with respect to the parameters. 

Preferably, the ophthalmic lens is a multifocal progressive lens. 

Preferably, the function includes a term involving the difference between 
the astigmatism induced by the ophthalmic lens and a predetermined astigmatism 
distribution. The predetermined astigmatism distribution describes different 
25 astigmatism than that required by a prescription for the ophthalmic lens. 

Preferably, the function includes a term involving a power other than 2 of 
the astigmatism induced by the ophthalmic lens. 

Preferably, the function includes a term involving the difference between 
the prism induced by the ophthalmic lens and a predetermined prism distribution. 
30 Preferably, the function includes a term involving the difference between 

the gradient of a characteristic and the gradient of a predetermined characteristic 
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distribution, where the characteristic is selected from a group including: power, 
astigmatism, and prism. 

Preferably, the function includes the square of the difference between the 
power induced by the ophthalmic lens and a predetermined power distribution and 
5 other terms related to the power induced by the ophthalmic lens. The other terms 
related to the power induced by the ophthalmic lens include terms of the form 

\p(l,m)-Py(l,mY for at least one value of p other than 2, where P(/,m) denotes 
the power induced by the optical element, P y (l,m) denotes the predetermined 
power distribution, and (l,m) parameterize the surfaces of the ophthalmic lens. 
io Preferably, the function includes a term related to the thickness of the 

ophthalmic lens. 

There is also provided in accordance with an embodiment of the present 
invention a method for designing at least one surface of a multifocal optical 
element. The method includes the step of concurrently considering in a design of 

15 the at least one surface both tone and progressive portions of a prescription for the 
multifocal optical element. 

There is also provided in accordance with an embodiment of the present 
invention a method for designing a surface of an optical element. The method 
includes the steps of prescribing initial wavefronts, selecting an initial 

20 representation of the surface, the representation including parameters, and 
precomputing eikonal functions between points in the vicinity of the initial 
representation and points in the vicinity of the initial wavefronts. When optimizing 
a cost function dependent upon the parameters, a refracted wavefront for each of 
the initial wavefronts is calculated from the precomputed eikonal functions. 

25 Moreover, the step of optimizing includes the step of analytically 

computing derivatives of the cost function using the precomputed eikonal 
functions. 

Furthermore, if during the step of optimizing, a current representation of 
the surface varies too much from the initial representation, then eikonal functions 
30 between points in the vicinity of the current representation and points in the vicinity 
of the initial wavefronts are computed. 
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There is provided in accordance with an embodiment of the present 
invention a method for designing a surface of an optical element. The method 
includes the steps of prescribing initial wavefronts, selecting an initial 
parameterized representation of the surface, choosing a function in the 
parameters, and optimizing the function with respect to the parameters. The step 
of optimizing includes the steps of calculating a refracted wavefront for each of the 
initial wavefronts and analytically computing derivatives of the function. 

There is provided in accordance with an embodiment of the present 
invention a method for designing at least one surface of an ophthalmic lens. The 
method includes the steps of representing the at least one surface by parameters, 
choosing an astigmatism distribution and an astigmatism direction distribution, 
choosing a function in the parameters and optimizing the function with respect to 
the parameters. The function includes a term involving the astigmatism 
distribution, the astigmatism direction distribution and the astigmatism induced by 
the ophthalmic lens. 

Preferably, the ophthalmic lens is a multifocal progressive lens. 

Preferably, the astigmatism distribution and the astigmatism direction 
distribution are determined from an eyeglass prescription. 

Preferably, the term is of the form 

%w l (!,my(ph([,m)-a 2l Q,m)Y -C*, (/,«)+ w t (f,m}a^(f,m)-C rA Q t m^ 

where (l,m) index the parameters for the at least one surface, C K>1 is the 
astigmatism distribution, C V1 is a second astigmatism distribution, a n , a n 
and a n are the coefficients in a quadratic expansion of a wavefront 
refracted by the lens in a coordinate system chosen for each (l,m) in 
accordance with the astigmatism direction distribution, w, and w 2 are 
weight distributions and p is an exponent. 
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BRIEF DESCRIPTION OF THE DRAWINGS AND APPENDIX 

The present invention will be understood and appreciated more fully from 
the following detailed description taken in conjunction with the appended drawings 
in which: 

5 Fig. 1 is a schematic illustration of an eye, and a lens having at least two 

surfaces, helpful in understanding the present invention; 

Fig. 2 is a flowchart illustration of a method for designing optical elements 
using wavefront calculations, according to an embodiment of the present invention; 
Fig. 3 is a schematic illustration of surfaces of a lens and of chief rays 
10 passing therethrough, according to an embodiment of the present invention; 

Fig. 4 is a schematic illustration of surfaces of a lens and of chief rays 
passing therethrough, according to an embodiment of the present invention, 
helpful in understanding the representation of localized wavefronts; 

Fig. 5 is a flowchart illustration of the optimization step of Fig. 2, according 
15 to an embodiment of the present invention; 

Fig. 6 is a schematic illustration of a chief ray and surfaces of an optical 
element, according to an embodiment of the present invention; 

Fig. 7 is a flowchart illustration of the eikonal method for an illustrative, 
non-limiting example, according to an embodiment of the present invention; 
20 Fig. 8 is a schematic illustration of the surface S? and a chief ray passing 

therethough, according to an embodiment of the present invention; and 

Appendix A is a derivation of equations for point objects and wavefronts 
used in the eikonal method of Fig. 7. 
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DETAILED DESCRIPTION OF THE PRESENT INVENTION 

The present invention is directed to a method for designing optical 
elements using wavefront calculations. An optical element is designed so that it 
converts a set of prescribed initial wavefronts into a set of associated wavefronts 
with desired properties. In general, it is not possible to solve the problem exactly, 
so an optimization method is used to produce the best optical element subject to 
the design constraints. The resulting designed surfaces of the optical element are 
integral surfaces and not composite surfaces. 

The present invention provides a method that enables design of a 
multi-surface optical element, in which any or all of the surfaces are designed, 
possibly simultaneously, using the method of the present invention. 

The present invention is directed to efficient optimization methods that 
provide the designer of the optical element with a large number of degrees of 
freedom to ensure a superior optical element. The degrees of freedom are the set 
of variables that define any of the surfaces of the optical element that need to be 
designed. It may be that some of the surfaces are fixed and cannot be changed. 
In that case, only the variables determining the non-fixed surfaces are considered 
to be degrees of freedom. In order to provide the designer with a large number of 
degrees of freedom, it is necessary to represent the surfaces efficiently and 
accurately. Applicants have found that it is advantageous to use several different 
surface representations. As will be described hereinbelow, some representations 
are useful for calculating wavefronts, while another representation is useful in an 
optimization process. 

In ideal situations, the surfaces are given explicitly in terms of well-known 
functions. Examples of such surfaces include spherical surfaces, aspherics, etc. 
In general, though, the designer should be allowed to design arbitrary surfaces, 
and therefore, in the present invention, discrete surfaces are used. Discrete 
surfaces are described in terms of finitely many points, for example, in the form of 
a set Q = {pt =(x / ,^,z / )}, / = 1, . . . , iST f where N is the total number of points. 

There are several efficient ways to describe discrete surfaces. Spline 
methods have the drawbacks mentioned in the background. U.S. Patent 
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Application Serial No. 09/262,341 to Kat2man et al. describes a finite element 
method for describing discrete surfaces. 

Applicants have discovered that a surface representation, termed 
"unconstrained patches", is highly efficient for the purpose of optical design. An 
arbitrary surface is described in terms of N discrete points, for example, in the 
form of the set Q mentioned hereinabove. In the vicinity of each discrete point 
Pi m (Wit*#) (" a Patch"). *e surface is approximated by a local function 
z = f t (x,y), where f, is a polynomial of a predetermined order in the variables x 
and y. The unknown polynomial coefficients are chosen to best fit the discrete 
values of the nearest neighbors to p, in the set Q. The number of neighbors for 
the fit is preferably taken to be larger than the number of unknown coefficients of 
/, to ensure accuracy and stability. This problem of determining the coefficients 
can be handled by techniques known in the art, such as the least squares method. 

Therefore, each point p, in the set Q is associated with a list of the 

coefficients of /,. It is noted that each patch is constructed separately and 

continuity of any kind between patches is not required. This method of 
"unconstrained patches" has several advantages. Firstly, the discrete set Q is 
arbitrary, and not constrained as in the finite elements method, or, even worse, as 
in the bicubic spline method. Secondly, the degrees of the polynomials /, are 

chosen by the designer, adding to the flexibility of the method. 

The present invention is particularly suitable for the design of ophthalmic 
spectacle lenses, including progressive addition lenses. Therefore, the detailed 
description that follows describes such an application, in which the initial 
wavefronts are refracted by the surfaces of the lens. However, it will be 
appreciated by persons skilled in the art that the present invention is equally 
applicable to the design of other optical elements including contact lenses and 
intraocular lenses and to the design of corneal shaping. 

Reference is now made to Fig. 1, which is a schematic illustration of an 
eye 100 and a lens 102 having at least two surfaces 104, according to an 
embodiment of the present invention. Eye 100 comprises a retina 106, an 
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eye-lens 108, a pupil 110 and a cornea 112. At least one of surfaces 104 of lens 
102 is unknown and will be optimally designed by the method of the present 
invention. 

An initial surface is chosen for each of the unknown surfaces 104 of lens 
102. Consider a ray 124A drawn from a central point in the eye towards lens 102 
at a spatial viewing angle of 6 . The lens designer specifies, for each viewing 
angle 9 , the distance that the wearer of the lens should be able to see in focus. In 
progressive lenses, for example, viewing angles above the horizon (shown in Fig. 
1 as a horizontal dotted line) are generally reserved for far vision, while viewing 
angles below the horizon are generally reserved for near vision. 

When ray 124A reaches lens 102, it is refracted by surfaces 104, whether 
they are known surfaces or initial guesses for the unknown surfaces. Therefore, 
the direction of ray 124A on the side of lens 102 away from the eye is generally 
different than the direction of ray 124A on the side of lens 102 near the eye. Ray 
124A is drawn until it reaches a sphere 115A centered at the central point in the 
eye and whose radius is the prescribed distance for the viewing angle 0 . Ray 
124A intersects sphere 1 1 5A at point 1 14A. 

If point object 114A were to emit, say, spherical wavefronts 120A, then 
they would be refracted by lens 102 into distorted wavefronts 122A. Since lens 
102 comprises initial guesses for the unknown surfaces, distorted wavefronts 122A 
are probably inconsistent with the wearer's prescription. In other words, the initial 
guesses for the unknown surfaces are probably incorrect An optimization process 
is used to determine the optimal surfaces 104 of lens 102, based on a cost 
function. During the optimization process, the unknown surfaces of lens 102 are 
varied from the initial surfaces, causing the direction of ray 124A on the side of 
lens 102 away from the eye to change, as indicated by arrow 126A. As the 
direction changes, the intersection of ray 124A with sphere 115A, namely point 
object 1 14A, changes. It will be appreciated, however, that since, by definition, 
point object 1 14A remains on sphere 1 15A, the lens designer's specification of the 
distance that the wearer of the lens should be able to see in focus is maintained 
throughout the optimization process. 
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The eye has a natural aperture - the eye's pupil 110. As the eye rotates 
while scanning the viewing field, only a small portion of each of the spherical 
wavefronts 120A emitted from the point objects actually gets transmitted through 
the pupil 110. This small portion of each spherical wavefront 120A is converted by 
lens 102 to a small portion of the corresponding refracted wavefront 122A. These 
small portions ("localized wavefronts") are centered around ray 124A, which is 
termed a "chief ray". Each point object 1 14 is associated with its own chief ray 124 
having a viewing angle 0 . 

Since the shape of refracted wavefront 122A contributes to the image 
seen by the eye, it characterizes the performance of lens 102 with respect to the 
point object 1 14A. Thus the designer wants the refracted wavefronts 122 to have 
desired properties in order to achieve high quality vision. The article "Testing and 
centering of lenses by means of a Hartmann test with four holes" by D. Malacara 
and Z. Malacara, published in Optical Engineering . 1992, vol. 31 no. 7, pp. 1551 - 
1555, explains how to read from the refracted wavef rant's shape the standard 
characteristics such as the power of the lens, its prism, its astigmatism, its coma, 
etc. A standard procedure in the ophthalmic industry is to define the vertex 
sphere. This is a virtual sphere 118 of a fixed radius known in the industry, for 
example 27 mm, centered at the center of the eye, having its vertex coinciding with 
the vertex of the back surface of lens 102. Commonly the wavefronts are 
propagated beyond the lens to virtual sphere 118, where the optical characteristics 
are computed. 

Reference is also made to Fig. 2, which is a flowchart illustration of a 
method for designing optical elements using wavefront calculations, according to 
an embodiment of the present invention. 

Initial surfaces are chosen for the unknown surfaces of the optical element 
(step 200). These initial surfaces are typically represented in one of the 
representations described hereinabove. Preferably, these initial surfaces are 
represented with a finite elements method or with the method of unconstrained 
patches. 

The designer of the lens then prescribes initial wavefronts (step 202). 
These initial wavefronts can be prescribed in any way. For simplicity, many point 
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objects 114 are selected, each emitting a spherical wavefront 120. In an 
embodiment, the designer specifies as a function of the spatial viewing angle 0 
the distances of the point objects from a central point in the eye. For example, in 
the design of a single vision lens, all the point objects are placed at a substantially 
fixed distance from the central point in the eye. 

The spherical wavefront 120 from each point object 114 propagates 
towards the eye. As it passes through lens 102, it is refracted by each of lens 
surfaces 104, thus resulting in a refracted wavefront 122. 

A chief ray is constructed, as explained hereinabove with respect to Fig. 1, 
for each of the point objects (step 204) using ray tracing techniques. 

The optimization process of the present invention requires repeated 
calculations of the intersection points of the chief rays with the surfaces of the 
optical element. When the surfaces are described discretely in a fixed global 
coordinate system, for example Cartesian coordinates or polar coordinates, as 
described hereinabove, this calculation is computationally intensive. Therefore, an 
alternative "ray representation", in which the surfaces are represented by their 
distance from predetermined fixed points along the chief rays, is preferred. The 
difference between the previous representation and the ray representation is 
analogous to the difference between the Eulerian and Lagrangian formulations that 
are well known in the theory of mechanics. 

A ray representation of the known and initial surfaces 104 is developed 
(step 206). Reference is now made briefly to Fig. 3, which is a schematic 
illustration of surfaces of the lens and of chief rays passing therethrough, according 
to an embodiment of the present invention. The lens consists of several surfaces 
S l9 S 2 ,.„ 9 S k . A chief ray can be identified by the viewing angle of the point object 
that it connects with the eye. The viewing angle is parameterized by discrete 
indices / and m. A chief ray r (l,m) t constructed as described hereinabove 
consists of several straight line segments whose directions change as the ray is 
refracted by the lens surfaces. The intersection point of chief ray r(l,m) with the 
surface Sj is denoted by g y (/,/w). The set {qj(l,m)} provides a discrete 
representation of the surface S, . This is the "ray representation". 
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Fixing the directions of the various segments of chief ray r(l,m) (using 
Snell's law), the intersection poMsq l (l,m),q 2 (l,m),...,q k (l,m) of chief ray r{l,m) 
with the lens surfaces can be represented by the sequence of distances between 
adjacent intersection points: f,(/,w)=|g, (l,m], t 2 {l,m)=\q 2 i},m)-q^,m\, .... 

t k (l,m) = \q k {l,m)-q k _ x {l,rn\. The distance / 1 (/,/ra)=|g,(/,m)| measures the 
distance along chief ray r(/,m) from starting point O (typically a central point in the 
eye) to the surface S x . Similarly, the intersection points of a chief ray r(l',m') also 
starting from O can be represented by the sequence of distances 

^>')=W>')|. '2('>0Hfc(/>0-<?.(/>% 

h(WhMW)-q k -,(l',m% 

The set of values {tj(l,m)} also provides a discrete representation of the 

surface Sj. The advantage of this Lagrangian description is that it automatically 

provides the intersection points of the chief rays with the lens surfaces. In the 
optimization process, the surfaces are modified in an iterative process that 
eventually converges to an optimal solution. With this Lagrangian description, the 
surfaces are modified by varying the values {*,(/, m)} along the chief rays. This 

has several advantages: firstly, there is no need to compute the intersection points 
with the chief rays, because these intersection points are now built in the surface 
representation; secondly, the surface is conveniently modified in only one 
dimension - along the chief ray. 

A ray representation can also be used in conjunction with other surface 
representations, for example in a finite elements method. Consider a system of 
rays parameterized by two discrete indices (l,m) and a single refractive surface 
S, . Define a reference surface S p also parameterized by (l,m) and set the 
distance between S p and 5, along any ray to be t(l,m) . 5, is fully determined by 
the fixed reference surface S p and by the function t. A finite element 
representation for the function t(l,m) is then constructed, as described in the 
classical texts of the finite element methods, and in the particular optical context by 
U.S. Patent Application Serial No. 09/262,341 to Katzman et al. 
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Once the surfaces have been represented by a ray representation, one 
needs to calculate the distortion of the initial localized wavefronts as they are 
refracted by the lens. A standard way to compute these distortions is to trace 
many rays through the lens. This method has the serious disadvantage that it 
requires the computation of the intersection of many rays with each of the lens 
surfaces, and the refraction of the rays at the intersection. Such calculations make 
this approach highly inefficient, since, in the optimization process, many 
computations of the distortion of localized wavefronts are made. 

Since, as mentioned hereinabove, only a small portion of each of the 
spherical wavefronts 120 emitted from the point objects 114 actually gets 
transmitted through the pupil 110, only the localized refracted wavefront 122, 
which is centered around the chief ray 124 leading from the point object 1 14 to the 
eye, needs to be calculated. The localized refracted wavefront can be used to 
compute the various wavefront characteristics by methods that are known in the 
art. 

The localized wavefronts are represented (step 208 of Fig. 2) using Taylor 
expansions. Reference is now made briefly to Fig. 4, which is a schematic 
illustration of surfaces 5 y of the lens and of chief rays r(l,m) passing 

therethrough, according to an embodiment of the present invention. A localized 
spherical wavefront 420 is refracted by the lens surfaces into a localized refracted 
wavefront 422 beyond the lens. It is necessary to calculate the localized refracted 
wavefronts. 

According to an embodiment of the present invention, the calculation is 
based on a direct computation of the first few terms of the Taylor expansion of the 
localized wavefront as it is refracted by an optical surface. For this purpose, 
"front-associated" coordinate systems are defined. These are local coordinate 
systems having a ^ axis pointing in the direction of the chief ray and variables £ 
and 7] in the orthogonal plane. Since the wavefront is orthogonal to the chief ray, 
there are no linear terms in the Taylor expansion for the wavefront's local 
representation in these front-associated coordinates. Thus the leading order terms 
in the Taylor expansion for the wavefront are given- by Equation 1 : 
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£ = a,,£ 2 + a, 2 #7 + <*vrf + • " ( Ec l- 1 ) 

where a,,, a, 2 and are Taylor coefficients. 

According to an embodiment of the present invention, only quadratic terms 
in the local Taylor expansion are considered. This effectively acts to restrict the 
wavefront to the portion that is transmitted through the pupil. These quadratic 
terms are related to the principal curvatures of the wavefront and their directions. 
The relations between the Taylor coefficients and the power and astigmatism of 
the wavefront are given by Equation 2: 

Power = a n + a n 

Astigmatism = -a n J +a? 2 

The variables £ and rj are defined above with respect to an arbitrary 
orthogonal coordinate system perpendicular to the chief ray. In some applications, 
as will be described in more detail hereinbelow, it is preferred to choose a specific 
coordinate system in the plane perpendicular to the chief ray. For example, 
Applicants have discovered that it might be preferable to rotate the ig,fj) 
coordinate system so that the new axes will be along the directions of the principal 
curvatures of the wavefront. It will be appreciated that in this rotated coordinate 
system, the Taylor coefficient a n , which contributes to the astigmatism as per 
Equation 2, has a value of zero. 

Therefore, the triplets { [a n {l.m\a n (},m\a n (l,m)]} represent the 

wavefront. The front-associated coordinate system thus leads to 

a representation 9 a{£° ,0%°] for the localized incoming wavefront, and the 
front-associated coordinate system (l^,^^,^^) leads to a representation 
[off* \an F \cc^\ for the localized refracted wavefront. It is these Taylor 
coefficients for the localized refracted wavefronts for each chief ray r(/,/w) that 
need to be calculated. 

Since the localized wavefronts are of finite extent, it is expected that they 
will overlap somewhat. A given area on the lens may focus wavefronts from one 
point object differently than it focuses wavefronts from another point object. 
Therefore it may happen that some lens shapes are advantageous for viewing 



15 



WO 01/92948 



PCT/IL01/00492 



some objects but not for viewing other objects. The optimization process 
compromises between the requirements of good lens performance for all the point 
objects. In order to compute a lens providing high quality vision, the designer 
determines a cost function to be optimized (step 210 of Fig. 2). 

In an embodiment of the present invention, the cost function G is written 
in the form given in Equation 3: 

M 

where (/,/») are the parameters used to describe the lens surfaces, D is the 
domain in which / and m are defined, w, are the weight functions used by the 
designer to emphasize different factors in the cost functions, S stands for the 
collection of lens surfaces S l} S 2 ,...,S k ,and the sum is over the M various factors 
U, . The weight functions w, can depend on the parameters / and m and on the 
surfaces S and their gradients. In an embodiment, some of the factors U t are 

chosen to be linear functions of the weights, in the form: 

U&m t S,w)=wfi,m t S)rft,m t S\ 

where T t are the factors in the cost function. 

The designer selects the factors U, and/or T, for optimization. In an 

embodiment the factors U, and/or T, are any combination of the following terms: 

• Required Power. For this purpose, the designer predetermines a function 
P r (f,m) that prescribes the desired power of the lens at the point (/,/n). 
Then the designer sets a term that is related to the deviation of the actual 
power P(l,m) from the desired power P Y (/,/«). In the prior art, the only 

power related term appearing in a cost function is |P(/,7«)-i^(/,»i)| 2 . 

Applicants have discovered that other power-related terms can improve the 
design of a lens. In particular, in the cost function, the power term could be 

a weighted combination of terms of the form \p(l,m)- P v {l,mf , for several 

choices of fi. 
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Required Astigmatism. Every progressive lens has some level of 
astigmatism C(l,m). In many applications, the designer attempts to 
minimize this astigmatism. For example, this is a design goal when 
designing semi-finished lenses or when designing a customized lens 
according to the specific prescription of a client where the prescription has 
no level of astigmatism. Such situations are termed "astigmatic-free" 
designs. In other applications, the designer wants to achieve a lens with a 
predetermined astigmatism distribution. For example, this is a design goal 
when designing a customized lens according to the specific prescription of 
a client where the prescription has some level of astigmatism. Such 
situations are termed "astigmatic" designs. It is appreciated that the 
astigmatism in astigmatic designs includes information on the required 
directions of astigmatism. 

Astigmatic-free designs can be accomplished, for example, by using 

a term of the form C(l,mf . Unfortunately, progressive lenses designed by 
minimizing the lens astigmatism in some parts of the lens often suffer from 
high disturbing astigmatism in other parts of the lens. Applicants have 
discovered that a lens design can be improved by other choices. In 

particular the form \C(l,mtf , for different choices of fi, is useful. One 
motivation for choosing a p other than 2, is that choosing, say, p = 1 gives 
extra weight to areas in the lens where the astigmatism is lower. In an 
embodiment, a weighted combination of terms of the form \C(l,m] fi for 
several choices of p is used in the cost function. 

In another preferred embodiment, a weighted combination of terms 
of the form \c(l,m)-C y (l,mY for several choices of p is used, where 
C v (},m) is an astigmatism distribution different from the client's prescription. 
In particular, even when the prescription has no astigmatism, a non-zero 
choice for C v (l,m) is recommended. It might seem odd initially to 

deliberately choose a non-zero C K (/,w), but Applicants have realized that 
by allowing low astigmatism at some parts of the lens, the designer has a 
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better chance of reducing the high, disturbing astigmatism at other parts of 
the lens. Furthermore, it is well known that the eye can adjust to low levels 
of astigmatism. 

Terms U l in the cost functions that relate to astigmatic designs will 
now be explained. The goal of an astigmatic design is to achieve 
predetermined distributions of astigmatism C v (l,m), and astigmatism 

directions p(l,m). Astigmatic designs are particularly useful in designing 
customized ophthalmic lenses according to a clients prescription. The 
prescription includes the prescribed power of the lens, the prescribed 
astigmatism and the direction of the prescribed astigmatism. In an 
embodiment, the prescribed astigmatism direction is used to select the 
distribution p(l,m) of astigmatism directions in the plane perpendicular to 
the chief rays. The required astigmatism values C y (l 9 m) are also selected 
according to the prescribed astigmatism. In an embodiment, the designer 
may select C v (l,m) and p(l,m) for each point (l,m) according to the 
prescribed astigmatism and prescribed direction regardless of the viewing 
angle. In another preferred embodiment, the distributions C y (l,m) and 

p(l,m) are determined from the prescribed astigmatism and the viewing 
angle by some rule. For example, a variation of Listing's law as described 
in the article by J. Porill, J.P. Ivins and J.P. Frisby, The variation of torsion 
with vergence and elevation", Vision Research, volume 39, pp. 3934-3950 
(1999) can be used for this purpose. 

In an embodiment, the astigmatism direction distribution is used to fix 

the (1,7) coordinate system in the orthogonal plane for each chief ray. 
Using this choice, the Taylor coefficient a ]2 in Equation 1 ought to have a 
value of zero, and therefore the astigmatism of the wavefront would be 
given by the difference between the Taylor coefficients ^ and . 

In another preferred embodiment, combinations of terms of the form 

w x (/, mU(a n (l 9 m) - a n (/, m)f - C 2 VtX (/, m)+ w 2 (/, m)a^ (/, m) - C va (l,mf 
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for various choices of p are used in the cost function, where w x (l y m) and 
w 2 (l 9 rn) are weights, the distribution C y l (l,m) is the desired astigmatism 

distribution according to the client's prescription and the distribution 
C K 2 (/,m) is selected by the designer to allow low levels of non-prescription 

astigmatism in various regions of the lens. Terms of this form have the 
advantage that during the optimization process, the designer can control 
both the required astigmatism direction (by minimizing the undesired 
coefficient a n ) and the astigmatism value (through the distributions 
Q,(/,m) and C K2 (/,m)). 

Cosmetic considerations. One of the criteria in the ophthalmic industry for 
the quality of a lens is its overall appearance, and in particular its thickness. 
Since the optimization methods of the present invention take the full 
structure of the lens into account, together with its constituent surfaces, 
Applicants have discovered that it can be useful to add to the cost function 
a term that is proportional to the lens thickness at the point (l,m). 

Required prism. The local behavior of each wavefront determines the local 
prism of the lens. In fact, the prism is already determined by the deflection 
of the chief ray. Applicants have realized that it is useful to predetermine a 
function Prism v (l 9 m) prescribing the desired prism for the lens at the point 

(/,m). In an embodiment, a weighted combination of terms of the form 

\PrismQ,m)- Prism y (l,my for several choices of p is used in the cost 

function. 

Variations in the lens characteristics. Applicants have realized that the lens 
performance can be improved by preventing too rapid variations in the lens 
characteristics. Alternatively, the designer may wish to obtain a prescribed 
distribution of the variations. For every lens characteristic, denoted by CH , 
the designer adds a weighted combination of terms of the form 

fiCH{Um)-VCH v (l,mY ', where V is the gradient operator, and 
VCH y (l 9 m) is the desired gradient distribution of CH . In an embodiment, 
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controlling the variations in the power, astigmatism and prism improves the 
lens quality. 

The design method based on the cost function G can be used to design 
surfaces that are progressive, astigmatic and prismatic. Moreover, any 
combination of the progression of the optical power, the astigmatism of the lens, 
and the prism of the lens are designed simultaneously and concurrently, thus 
leading to integral surfaces. It will be appreciated that the design of integral 
surfaces is superior to the prior art, as described in U.S. Patent No. 2,878,721 to 
Kanolt and U.S. Patent No. 6,019,470 to Mukaiyama et al., where two separate 
surfaces are designed, one for the optical power progression, and one for the 
required astigmatism, and then composed to obtain a composite surface. 
Consequently, in a lens with a composite surface, one or more of the constituents 
of the composite surface may damage some of the optimal properties of the other 
constituents. While in the design of a composite surface there is no correlation 
between the two parts (progressive and toric), the design of an integral surface 
takes both of them into account simultaneously, thus ensuring an overall 
optimization. In another example, the thickness of a toric contact lens is increased 
in certain regions of the lens so that when worn, it will stay on the eye in a 
particular orientation. The prism of the contact lens needs to be adjusted in those 
regions in order to compensate for the added thickness. Thus, it is beneficial to 
design one or more integral surfaces using a design method that takes the 
astigmatism and the prism into account simultaneously and concurrently. In a 
further example, there are design requirements for lenses having progressive 
power, astigmatism and prism. The design method described herein enables the 
design of an integral surface while concurrently taking into account all of these 
requirements. 

Once the cost function has been chosen, it needs to be optimized (step 
212 of Fig. 2). The optimization process consists of searching for those surfaces 
that yield the best value for the cost function. The cost function is optimized 
subject to the requirement that the surfaces never intersect and possibly subject to 
additional constraints given in the form of prescribed values for the surfaces and/or 
their derivatives at a selected set of points. The result of the optimization process 
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is a numerical representation of all designed surfaces of the optical element (step 
214). It will be appreciated that these designed surfaces are integral surfaces and 
not composite surfaces. 

Many optimization methods are known in the art. For example, common 
techniques include the conjugate gradient method and the Newton method. 
Essentially all methods are iterative. Reference is now made to Fig. 5, which is a 
flowchart illustration of the optimization step 212 of Fig. 2, according to an 
embodiment of the present invention. 

When the Rih step in the iteration begins, each unknown surface Sj has 

the shape S* (step 500). Localized spherical wavefronts emitted by point objects 

are refracted by the surfaces of the lens into localized refracted wavefronts beyond 
the lens. It is necessary to calculate the localized refracted wavefronts. 

In an embodiment, the Taylor coefficients for the localized refracted 
wavefronts beyond the optical element are calculated using a localized front 
refraction method (step 502). Referring back to Fig. 4, in this method, another 

local coordinate system (gf^tff* is constructed at the intersection point of 

the chief ray with the surface S k (farthest from the eye), so that its £ axis points in 

the direction of the normal of the surface. The surface S k is locally represented by 

its Taylor expansion in this coordinate system [a™ >a™ 9 a™]. By differentiating 

Snell's law, the Taylor coefficients of the localized wavefront after refraction by the 
surface S k can be expressed in terms of the Taylor coefficients of the incoming 

localized wavefront and the Taylor coefficients of the surface at the intersection 
point This localized front refraction method is repeated for each surface until the 
Taylor coefficients of the localized refracted wavefront beyond the optical element 
have been calculated. 

A formulation of the localized front refraction method is given in J. Kneisly, 
"Local curvature of wavefronts in an optical system", JOSA 54, 1964, pp. 229-235. 
A formula is developed that directly connects the principal curvatures and 
directions of a refracted wavefront with the principal curvatures and directions of 
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the wavefront before refraction and the principal curvatures and directions of the 
surface at the point of refraction. 

Applicants have discovered a novel method, based on eikonai functions, 
for calculating the localized refracted wavefronts beyond the optical element. 
Eikonai functions are well known classical objects in optics. O.N. Stavroudis, The 
Optics of Ravs. Wavefronts and Caustics . Academic Press, 1972, pp. 200 - 207 is 
a general reference on the history of eikonals and their applications in optics. The 
novel eikonai method involves precomputing eikonai functions and geometrical 
quantities related to the known and initial surfaces. It provides an efficient explicit 
computation of the localized refracted wavefronts, and provides an explicit 
computation of the coefficients [a/f '.a^ ,a^\ in terms of the free parameters. 
The eikonai method (step 504) is described in detail hereinbelow. 

The calculation of the shape of the localized wavefront after refraction by 
all surfaces of the optical element, using the eikonai method or the localized front 
refraction method or some other suitable method, results in localized refracted 
wavefronts beyond the optical element (step 506). 

Once the localized refracted wavefronts have been calculated, the cost 
function can be evaluated (step 508) and checked (step 510). If the cost function 
is sufficiently optimized according to some predetermined criteria, then the 
surfaces S* are optimal solutions for the unknown surfaces Sj (step 512). If, 

however, the cost function is not sufficiently optimized, then the surfaces Sf are 

updated for the next iteration (step 514) and the method continues the next step in 
the iteration (step 500). 

Each iteration step requires the calculation of the derivative of the cost 
function with respect to all its variables. In fact, popular optimization methods, 
such as the Newton method, require also the computation of the second 
derivatives. Since a flexible design includes a very large number of free variables, 
and since the computation of derivatives is required for each iteration, it is 
desirable to look for efficient methods of computing these derivatives. 

One possibility is to compute the derivatives numerically. For this 
purpose, one slightly varies a certain parameter and evaluates the resulting 



22 



WO 01/92948 



PCT/IL01/00492 



change in the cost function. Since the lens surfaces and the refraction of the 
wavefronts by them are characterized by local properties, there is no need to 
recompute the transfer of all of the wavefronts for each such small local variation. 
This localization accelerates the computation of the derivatives, and hence the 
entire optimization process. 

Another possibility is to compute the derivatives analytically, which can 
accelerate the optimization process dramatically. Applicants discovered that the 
eikonal method (step 504) enables this analytical computation of derivatives, as 
will now be described in detail using an illustrative, non-limiting example. The cost 
function involves power, astigmatism, prism, and thickness of the lens and the 
gradients thereof. We have already shown that power and astigmatism can be 
expressed in terms of the Taylor coefficients as in Equation 2. The prism can also 
be expressed in terms of the Taylor coefficients. The eikonal method enables the 
Taylor coefficients to be expressed explicitly in terms of the problem parameters. 
Clearly the thickness of the lens is also expressed explicitly in terms of the problem 
parameters. Therefore the derivatives of the cost function can by analytically 
computed. Reference is now made to Fig. 6, which is a schematic illustration of a 
chief ray r(l,m) and surfaces S u S 2 ,... 9 S k of an optical element, according to an 
embodiment of the present invention. In this specific example, the surfaces 
S 2 ,...,S k are assumed to be known and fixed, while the surface S { is to be 
optimized with respect to a cost function. The refraction index in the region 
between S x and the eye is denoted « 0 , and the refraction index in the region 

between S x and S 2 is denoted n x . Reference is also made to Fig. 7, which is a 
flowchart illustration of the eikonal method for the current non-limiting example, 
according to an embodiment of the present invention. 

' The design method begins with an initial choice for Si denoted by S* 
(step 700). Each chief ray r(l,m) is traced from an origin O (a central point in the 
eye) through the optical element until a point object o 0 (/,/w). The point object 
o°(l 9 m) lies on a sphere 600 centered at the origin O and having a radius of fixed 
distance \o{l 9 m\. This distance is fixed by the viewing angle of the chief ray 
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r(l,m). The distance from the origin O to the intersection point of the chief ray 
r(l 9 m) with iSf is denoted t*(l 9 m). The direction vector of the chief ray r(l,m) is 
denoted £(/,w). The set {t°(hm\i(l,m)} fully provides a discrete form of the initial 
surface S? . The choice of initial surface S° fixes the initial distances t*(l,m) and 
the initial point objects o°(l,m). 

For each chief ray r(/, m) a local orthogonal coordinate system having one 
of its axes in the direction of is defined (step 702). It will be appreciated 
that this is the local coordinate system if*" of Fig. 4. It should be 

noted that this coordinate system is independent of the value of t*(l,m), and 
therefore it will remain fixed throughout the optimization process. 

The eikonal function E is classically defined as the optical distance 
between two points. For the spherical wavefront emitted by the point object 
o°(l, m), the shape of the refracted wavefront at the surface 5° can be calculated 
from the eikonal function between o°(l,m) and t°(l+p,m + q) for all combinations 
of ^,# = 0,11. Thus in this preferred embodiment, each point object o°(l,m) has 
nine eikonal functions, (/ + p,m + q\o°(l 9 m)) for all combinations of p,q = 0,±1 , 
associated with it Since the surfaces S 29 ... 9 S k are known, Fermat's principle (also 
known as the principle of least time) can be used to calculate the distance of the 
optimal optical path (i.e. the eikonal function) between o°(l,m) and tf(l+p,m + q) 
for all combinations of p,q = 0,±1 (step 704). 

In the optimization process, the surface S x might deform from its initial 
shape 5,°. In other words, the initial values t*(l,m) might change into, say, 
t x (},m). Consequently, the associated point objects o°(l,m) might vary too into 
o(l 9 m). However, since the chief ray r(l,m) below the surface S x has not changed 
its viewing angle, the new point object is confined to the sphere 600 of radius 
|o(/,m). 

It is necessary at each step of the optimization process to calculate the 
shape of the refracted wavefront at the surface 5, due to the spherical wavefront 
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emitted by the point object o(l,m). This will be accomplished by calculating the 
eikonal function between o(l,m) and t x (l+p,m + q) for all combinations of 
p,q = 0,±1 (nine eikonal functions in all). If surface S x is not too deformed from the 
initial surface Sf , then the eikonal function between o(l,m) and a particular 
t x {l+p,m + q) can be expressed using a local Taylor expansion in the following 
form: 

E(tr,o)= E(t?(l + P> m + q\o°(l,m)) 

+iz! l 7 p - 9 (°j- o0 A l > m )) 

+ % ^™( t ™-t Q l (l + p i m + qtf . (Eq.4) 

+ tziS**(pi -o%m)loj ~o%m)) 

where the Taylor coefficients are denoted by ziT**' ^iT™. ^j^™. 

ifr". xfr** z l g p "> z l sT q > z& M > z%?**x\%™> &« and *™ denotes 

t x (l+p,m + q), o denotes o(l,m) and o } denotes the j-th component of o . 

However, the Taylor coefficients need to be determined in advance. 
Fermafs principle is used to calculate the value of the eikonal function between 
any point o in the vicinity of o°(l 9 m) and any point tf* in the vicinity of 
t*(l+p,m + q) (step 704). For each combination of /?,# = 0 9 ±1, these eikonal 
values are then substituted into Equation 4, thus providing systems of equations 
for the unknown Taylor coefficients that can easily be solved, for example by the 
least squares method. The number of eikonal values calculated in step 704 is 
greater than the number of unknown coefficients, to ensure accuracy and stability. 

The design process consists, therefore, of finding for each (l,m) the 
values ti(l,m) along the chief rays that optimize the cost function and the 
corresponding point objects o(l 9 m). Eventually it will be shown that the point 
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object o(l 9 m) is determined by the t x (l + p,m + q) values. So, in fact, the free 
variables for optimization are the points in the set t x (l,m). 

The optimization process is iterative. Consider the Rth step in the 
iteration. The surface S x has the shape S R parameterized by the collection 

t R (l 9 m) (step 706). For each l,m, the point object o R (l 9 m) associated with the 
chief ray whose direction vector is r (l,m) needs to be computed (step 708). The 
computation of step 708 uses the local coordinate system defined in step 702. 
Having found o R (l,m), Equation 4 is used to compute £(/,*(/ + p,rn + q\o R {},m$) 
for all combinations of /?,g = 0,±l (step 710), and from these values, the local 
shape of the wavefront is computed (step 712), which in turn is used to evaluate 
the cost function (step 714). If the cost function is not sufficiently optimized (step 
716) then the method continues from step 706, with the next iteration of the 
optimization process. If, however, the cost function is sufficiently optimized (step 
716) according to some predetermined criteria, then the shape S R parameterized 
by the collection t R (l,m) defines the optimized surface S x (step 718). 

Step 708, in which the new point object o R (l,m) for each chief ray is 
computed, will now be explained in more detail. The new point object 
could be found by tracing the chief ray r (/,/«) through the modified optical 
element. But this will result in a very complex dependence of o R (l,m) on the free 
variables t R (l 9 m). Therefore, an alternative approach is used which utilizes the 
local nature of the entire problem and, in particular, the precomputed eikonal 
function, to develop the following simple equation for the o R (l,m): 

e^^-n^'", (Eq.5) 
where the 2-dimensional vectors s Lm and <x> Km are defined as follows: 

^= ^(f+P 9 m + q yie(P.q)^ (Eq-6) 

p t q=±l t 0 

and 

M-±1.0 

for n m 1,2 . 
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The vectors e lJn and <o ljn involve the following fixed geometrical 
quantities: 

tfHp.9) = fQ.my("(t + P.m + q) (Eq.8) 
*(lW)-£~f + A« + *K'M 
and also involve the quantities wf* and ef£ , which are defined as follows: 

and 

C = 4* ( / + P>m + q\o%m))-E(t*Q 9 ml^ (Eq. 10) 

Equation 5 is an equation for the new point object o R (l,m) through the 

eikonal functions of Equation 10. 

The coefficients of the local wavefront are given explicitly by: 

a lm = [a'' w \ (n 0 + E l " ) , (Eq. 1 1 ) 

where a 1 * is the three-dimensional vector whose entries are the Taylor 
coefficients fan 9 a£ 9 a*£) of the refracted wavefront. Equation 11 involves the 
following 3x3 matrix A /tOT and three-dimensional vectors W /,m and E /,m : 



(^J^m J 2 / 4 2] A^^A^^t 4 ^ Alf tP,q A^^ q t 4 

W=0.±1 P*$=0,±1 P,9«0.±l 



_p,g=0,±l Ptf=0.±l M=0.±1 

where f 4 is a shortened notation for (if (/ + p 9 m + q)J , 



, (Eq.12A) 



p,?=0,±l 



(Eq. 12B) 
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(Eq. 12C) 



where the canonical quantities 4?' M >4$' P *>4"' P ' 9 are defined as follows: 



Equation 5 and Equation 11, which are derived in Appendix A, are the 
fundamental equations in the novel eikonal method. They provide explicit 
expressions for the point objects o R (l,m) and for the form of the local wavefront for 
each chief ray r(l,m). The equations contain a number of geometrical quantities 
such as 4g,6g.*S and A^- M ,A^' M ,4a'"' 9 ^at can be precomputed at the 

initial stage and then stored in memory. 

It might happen that during the iterative optimization process the surface 
5, will deform appreciably from its initial shape S?, and therefore the use of 
Equation 4 for the calculation of the local wavefront is no longer appropriate. 
Applicants have discovered that, in this case, it is preferable to recompute the 
eikonal function for points in the vicinity of the new surface 5, and to restart the 
process. It should be noted, though, that the eikonal function is computed only a 
small number of times in the iterative process. 

It will be appreciated by persons skilled in the art that the present invention 
is not limited by what has been particularly shown and described herein above. 
Rather the scope of the invention is defined by the claims that follow Appendix A: 



4r- 9 =te(/>>*)) 2 

A£">< =(b$(p,qjf 



(Eq.13) 
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APPENDIX A 

This appendix is a derivation of Equations 5 and 11, relating the point 
objects o R (l,m) and the coefficients of the refracted wavefronts to the free 
variables t R (l, m) for the J?th step in the iteration of the optimization process. 

Reference is now made to Fig. 8, which is a schematic illustration of the 
surface S R and a chief ray r(l,m) passing therethough, according to a preferred 
embodiment of the present invention. The chief ray intersects a sphere 800 at the 
point object o R (l,m). 

In order that the wavefront originating at the point object o R (l,m) 
propagates in the direction of the chief ray r(/,m), and thus is exactly aimed at a 
central point O in the eye, it must be that the normal to the wavefront at *,*(/, /w) is 
in the direction of f(/,/n). This condition is translated to a condition on the point 
object o R (l,m). For this purpose the plane PL orthogonal to f(l,m) at t R (l,m) is 
considered. The eikonal functions for points on that plane are evalulated. 

The eikonal function E in the plane PL is now expressed in the 
orthogonal local coordinate system ) that is already defined there. 

The restriction of E to the plane PL in these coordinates is denoted by a(^,rj). 
cr(t;,T}) is evaluated by solving the classical eikonal differential equation. The 
values of E at the points + p,m + q) for p,q = 0,±1 are the initial conditions for 
the differential equation. Equivalent^, <r(£ t tj) can be computed by adding to the 
above-mentioned initial values the optical distance along the rays emanating from 
each of the t R (l + p,m+q) points until the rays intersect the plane PL. 

A simple geometrical calculation then provides that the orthogonal 
projection of the points t*(l + p,m + q) onto the plane PL is given by: 

[$%,Ti?2)=t R Q +P,m + qXb?i<j>,q),b< 1 J<j> i qy) , (Eq. A1) 
where the geometrical quantities bfLb^bfj, are defined in Equation 8: 
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b%(p,q) = i* EF (l,rnyC* BF (l + P,rn + q) 

The condition on the point object o R (l,m) is that the normal to the 
wavefront at t R (l,m) is in the direction ^(^m). This is equivalent to the 
condition: 



da_ 
H 



(0.0) 



da 
drj 



= 0. (Eq. A2) 

(0,0) 



The Euclidean distance between a point t?(l + p,m + q) and its projection 
onto the plane PL is given by: 

"tf =t?(l + P,m + q)bS(p,q)-t?(l,m) (Eq. 9) 

It follows that, to linear order, the optical distance between point 
t?(l + p,m + q) and the plane PL is n 0 wf*. Since the condition of Equation A2 

holds, to linear order, the following equation is obtained: 

<rtftf,rif*)-cr(0,0) = 0. (Eq. A3) 

In particular, the equivalence holds for the projected points defined in 
Equation A1: 

V%)-Eit?{lm\o%m))=Q . (Eq. A4) 

The left side of Equation A4 is precisely the difference between the eikonal 
values at t?(l + p,m + q) and at t?(l,m), plus the optical distance between 
t?Q + P> m + q) ar, d its projection onto the plane PL. Thus using the optical 
distance calculated above, the following equation is obtained: 

+ -0. (Eq.A5) 

where wf* is defined in Equation 9, and ef£ is defined in Equation 10: 

etf=E(t l R (l + p,m + q\o%m))-E{t l %m\o%m)). (Eq. 10) 

Equation A5 is a system of 9 equations for the unknown point object 

o R (l,m). There are many ways to stabilize it, for example, a least squares 
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approximation. An alternative, natural way to solve the overdetermined system of 
Equation A5 is to replace Equation A3 by the following: 



da 



do- 



(Eq. A6) 



(0.0) 



(o.o) dT f 

Equation A6 is applied at the projected points, and using the linear-order 
optical distance, the following is obtained: 

d A + 

da 



(f?(l + p,m + q)bV(p,qj)- 



(0.0) 



(Eq. A7) 



d 1(.o,o) 



where p,q = 0,±1. 

Eq. A7 forms an overdetermined system of 9 equations in the two 



, do 
unknowns — 

a* 



(0.0) 



, 5<r 
and — 



(considering for the moment the unknown point 



(0.0) 



object o R (l,m) as a parameter). This overdetermined system can be converted 

naturally to a 2x2 system by the method of least squares. 

The normal equations for the least squares method can be written 
explicitly. It is convenient at this stage to introduce for each quadruple l,m,p,q 

the canonical 2x2 matrices: 

BSKlft" = b$(p,q)b}£<J>,q) , (Eq. A8) 

for i, j = 1,2. 

The matrix B'- m is also defined. Its ij entry is given by: 

B'J= ^(l + P) m + q)jBSKU^. (Eq. A9) 

The normal equations for the least squares approximation of Equation A7 



are: 



where 



5 / ' m Vo'(0,0) = « 0 ^+^ w f 



(Eq. A10) 



(0.0) 
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and where the 2-dimensional vectors s ,jn and a)'" are defined in Equations 6 and 
7, respectively: 

I l t?(l+P>m + q>blS(P><iye% (Eq. 6) 

p, ? =±l,0 

and 

for 77 = 1,2. 

Imposing now the condition Vo*(0,0) = 0 , Equation 5 is obtained for the 

two components of the point object o R (l,m): 

e^^-rw)**. (Eq.5) 

It is useful at this stage to point out certain computational aspects. Firstly, 

the geometrical matrices BSKL^*™ and B ljn are not used at all. Secondly, the 

right side of Equation 5 is a quadratic explicit expression in the variables t R , and 
does not depend at all on the point objects. Thirdly, the geometric quantities 
b Z> b i 2 J> b t 3 J depend only on the directions of the chief rays and are therefore 
computed only once (at the initial stage) and stored in memory. Finally, the 
vectors s ljn depend on the point objects o R {l,m) quadratically through the local 
expansion of the eikonal function. Therefore, the point object o R (l,m) can be 
expressed explicitly in terms of the variables t R (l 9 m) . 

Determining the point object o R (l,m) via Equation 5 implies that the linear 
terms in the expansion of the eikonal cr(£, 77) vanish. Therefore, to leading order, 
the local expansion takes the form: 

o*(£^)-cr(0,0) = £^ . (Eq.A11) 

The local quadratic expansion of Equation A11 is precisely the local expansion for 
the wavefront passing through t R (l 9 m). 

By substituting Equation A1 into Equation A11 and using the fact that 
<r(4^,Tj^)-cr(0 y 0)^n 0 w^ , one obtains the following: 
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W% + e% + P>>" + <l)}- 



+ bj»(p,q)b%(j>,q)ct l2 (l,m) .(Eq.A12) 

By defining the following canonical quantities as in Equation 13 in the detailed 
description: 

Equation A12 can now be written compactly in the following form: 

nXZ+e™=(t?{l + p,m + q))\4f< p «a\?+4? P M^ +A^a^],(Eq. A13) 
for p,q = 0,±l. 

Eq. A13 forms an overdetermined system of 9 equations in 3 unknowns 
(a,j (/,/»)). In order to write down the normal equations for the least squares 

approximation for its solution, the following 3x3 matrix a'-" 1 and three-dimensional 
vectors W'-" 1 and E ,,m are defined in Equation 12 of the detailed description: 



,g=04:l p>q*=Q>±l /?,$=0,±1 



p % q=QjA 



p,q=Q,±l 



_p,q=0,±\ m=o,±i m=o,±i 

where f 4 is a shortened notation for (/*(/ + p 9 m + , 

I i {tn^p^g)yA^w^ 



(Eq. 12A) 



p,9=0,±l 



^(rf(/ + ^/« + g)) 2 4; 

p,g=0,±l 



hm*p v q w P4 



p>i w p><! 

rr l,m 



and 



(Eq. 12B) 



33 



WO 01/92948 



PCT/IL01/00492 



E' ,m = 



Ebtft + P.m + qjfA}?"*, 



PA 



P£ 0 P4 



p,q=O t ±\ 



P4 0 P4 



(Eq. 12C) 



The normal equations are then given by: 

A^a'" = « 0 W + E (Eq. A14) 

where a l/tt is the three-dimensional vector whose entries are the Taylor 
coefficients ,or,'f of the refracted wavefront. The coefficients of the local 
wavefront are given explicitly by Equation 1 1 of the detailed description: 

a i, m = [ A /.m Y ( floW + E ) t (Eq. 11) 
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What is claimed is: 

1 . A method for designing at least one surface of an optical element, 
the method comprising the steps of: 

for each particular surface of said at least one surface: 

creating a representation of said particular surface with a set of 
discrete points; and 

for each discrete point in said set: 

selecting nearest neighbor points from said set; 
approximating said particular surface in a vicinity of said discrete 
point by a polynomial of a predetermined order in two variables of 
said particular surface; and 

determining coefficients of said polynomial according to said 
nearest neighbor points without requiring continuity between 
polynomials for neighboring discrete points. 

2. A method for designing at least one surface of an optical element, 
the method comprising the steps of: 

representing said at least one surface using the method of unconstrained 
patches; 

choosing a function in parameters of said at least one surface; and 
optimizing said function with respect to said parameters. 

3. A method for designing at least one surface of an ophthalmic lens, 
the method comprising the steps of: 

representing said at least one surface with parameters; 
choosing a function in said parameters, said function including a term 
involving the difference between the astigmatism induced by said 
ophthalmic lens and a predetermined astigmatism distribution, said 
predetermined astigmatism distribution describing different astigmatism 
than that required by a prescription for said ophthalmic lens; and 
optimizing said function with respect to said parameters. 

4. The method of claim 3, wherein said ophthalmic lens is a multifocal 
progressive lens. 
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5. A method for designing at least one surface of an ophthalmic lens, 
the method comprising the steps of: 

representing said at least one surface with parameters; 

choosing a function in said parameters, said function including a term 

involving a power of the astigmatism induced by said ophthalmic lens, 

said power having a value other than 2; and 

optimizing said function with respect to said parameters. 

6. The method of claim 5, wherein said ophthalmic lens is a multifocal 
progressive lens. 

7. A method for designing at least one surface of an ophthalmic lens, 
the method comprising the steps of: 

representing said at least one surface with parameters; 

choosing a function in said parameters, said function including a term 

involving the difference between the prism induced by said ophthalmic 

lens and a predetermined prism distribution; and 

optimizing said function with respect to said parameters. 

8. The method of claim 7, wherein said ophthalmic lens is a multifocal 
progressive lens. 

9. A method for designing at least one surface of an ophthalmic lens, 
the method comprising the steps of: 

representing said at least one surface with parameters; 
choosing a function in said parameters, said function including a term 
Involving the difference between the gradient of a characteristic and the 
gradient of a predetermined characteristic distribution, said characteristic 
selected from a group including: power, astigmatism, and prism; 
optimizing said function with respect to said parameters. 

10. The method of claim 9, wherein said ophthalmic lens is a multifocal 
progressive lens. 

11. A method for designing at least one surface of an ophthalmic lens, 
the method comprising the steps of: 

representing said at least one surface with parameters; 
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choosing a function in said parameters, said function including the square 
of the difference between the power induced by said ophthalmic lens and 
a predetermined power distribution and other terms related to the power 
induced by said ophthalmic lens; and 
optimizing said function with respect to said parameters. 

12. The method of claim 11, wherein said ophthalmic lens is a multifocal 
progressive lens. 

13. The method of claim 11, wherein said other terms related to the 
power induced by said ophthalmic lens include terms of the form 

\P(l,m) -Py(l 9 mY for at least one value of p other than 2, where P(l 9 m) 

denotes the power induced by the optical element, P v (l 9 m) denotes the 

predetermined power distribution, and (l 9 m) parameterize the surfaces 

of said ophthalmic lens. 

14. A method for designing at least one surface of an ophthalmic lens, 
the method comprising the steps of: 

representing said at least one surface with parameters; 

choosing a function in said parameters, said function including a term 

related to the thickness of said ophthalmic lens, and 

optimizing said function with respect to said parameters. 

15. The method of claim 14, wherein said ophthalmic lens is a multifocal 
progressive lens. 

16. A method for designing at least one surface of a multifocal optical 
element, the method comprising the step of: 

concurrently considering in a design of said at least one surface any 
combination of an astigmatism distribution, a progressive power distribution 
and a prism distribution. 

17. A method for designing a surface of an optical element, the method 
comprising the steps of: 

prescribing initial wavefronts; 

selecting an initial representation of said surface, said representation 
including parameters; 
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precomputing eikonal functions between points in the vicinity of said initial 
representation and points in the vicinity of said initial wavefronts; and 
when optimizing a cost function dependent upon said parameters, 
calculating from said precomputed eikonal functions a refracted wavefront 
for each of said initial wavefronts. 

18. The method of claim 17, wherein said step of optimizing comprises 
the step of: 

analytically computing derivatives of said cost function using said 
precomputed eikonal functions. 

1 9. The method of claim 17 further comprising the step of: 

if during said step of optimizing, a current representation of said surface 
varies too much from said initial representation, computing eikonal functions 
between points in the vicinity of said current representation and points in the 
vicinity of said initial wavefronts. 

20. A method for designing a surface of an optical element, the method 
comprising the steps of: 

prescribing initial wavefronts; 

selecting an initial representation of said surface, said representation 
including parameters; 

choosing a function in said parameters; and 

optimizing said function with respect to said parameters, said step of 
optimizing including the steps of: 

calculating a refracted wavefront for each of said initial wavefronts; and 
analytically computing derivatives of said function. 

21. A method for designing at least one surface of an ophthalmic lens, 
the method comprising the steps of: 

representing said at least one surface by parameters; 

choosing an astigmatism distribution and an astigmatism direction 

distribution; 

choosing a function in said parameters, said function including a term 
involving said astigmatism distribution, said astigmatism direction 
distribution and the astigmatism induced by said ophthalmic lens; and 
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optimizing said function with respect to said parameters. 

22. The method of claim 21, wherein said ophthalmic lens is a multifocal 
progressive lens. 

23. The method of claim 21, wherein said astigmatism distribution and 
said astigmatism direction distribution are determined from a 
prescription for an eyeglass. 

24. The method of claim 21 , wherein said term is of the form 

^Amy(a n (hm)-a a (l,m)f -C* I (/,»i)+ w 2 (/,m)4(/,m) - C rA (l,m^ wh 

ere (l,m) index said parameters for said at least one surface, Q , is 
said astigmatism distribution, C V2 is a second astigmatism distribution, 
a n , a n and a n are the coefficients in a quadratic expansion of a 
wavefront refracted by said lens in a coordinate system chosen for each 
(l,m) in accordance with said astigmatism direction distribution, w, and 
w 2 are weight distributions and p is an exponent. 
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